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Abstract 



We study the distribution of the time to explosion for one-dimensional diffusions. We relate this ques- 
tion to computing the expectations of suitable nonnegative local martingales, and to the distributions 
of related diffusions with unit variance. Moreover, we characterize the distribution function of the 
time to explosion as the minimal solution to a certain Cauchy problem for an appropriate parabolic 
differential equation; this leads to alternative characterizations of Feller's criterion for explosions. 
We discuss in detail several examples for which it is possible to obtain analytic expressions for the 
corresponding distribution of the time to explosion, using the methodologies developed in the paper. 
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1 Introduction and summary 

X ■ 

j_j '. Precise conditions for whether or not a one-dimensional diffusion process explodes in finite time have 

been developed, most notably by William Feller. To the best of our knowledge, and rather surprisingly, 
the distribution of the explosion time has rarely - if at all - been the subject of investigation. With this 
work we hope to help close this gap or, at the very least, to narrow it. 

We start by introducing the framework of this pape r and by discuss ing the Feller test in Section [2] 



Then, in Section [3j we recall and generalize a result of iMcKeanl (119690 that associates the distribution 



of an explosion time to the expectation of a related nonnegative local martingale. This approach can be 
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considered a tr ansformation of the probability me asure. We a l so dis cuss a transformation of space, based 
on the ideas of lLampertil (119641) . iDossI (119771) and lSussmann dl978l) . 

We study analytic properties of the distribution function in Section |4j and derive alternative char- 
acterizations of Feller's criterion for explosions. In particular, we characterize the tail P^(5 > t) of 
the probability distribution function of the time-to-explosion S, viewed as a function of time t and 
starting position £, as the smallest nonnegative solution of an appropriate partial differential equation of 
parabolic type. Finally, in Section [5] we provide several examples, to which we apply the methodologies 
of Sections Sand H 

Although there does not seem to be much work on the distribution of the time to exp losion, there ex- 
ists a h uge literature on the computation of first-passage times by diffusions. We refer to Pitman and Yor 
(|2003h and Section 2 in lGoing-Jaeschke and Yorl(|2003|) for several pointers to this literature. 



2 Some essentials of one-dimensional diffusions 



?, r) with — oo<£<r<oo and consider the stochastic differential 



We fix an open interval / = 
equation 

dX(t) = s(X(t))(dW(t) + b(X(t))dt) , X(0) = £, 

where £ G I and W(-) denotes a Brownian motion. We shall assume that b : / 
M \ {0} are given measurable functions that satisfy the condition 



and s 



K 



1 



s 2 (y) 



*(v) 



dy < oo, for every compact set K C I. 
In other words, we are assuming that both l/s 2 (-) and the "local mean/variance ratio" function 

f(-: 



(2.1) 

I -> 

(2.2) 



KO _ b(-M-) 
*(■) 



*»(■) 



(2.3) 



are locally integra ble over I. From the same arguments as in lEngelbert and Schmidt! (|1984i Il99lh or 
Theorem 5.5.15 in lKaratzas and Shrevel (1 199 lb . it follows that the stochastic differential equation (12.11 ) 
admits a weak solution, unique in the sense of the probability distribution and defined up until the "ex- 
plosion time" 



S := inf {t > : X(t) <£ (I, r)} = lim S n 

n\oo 



S n :=inf{i>0:X(t) ^ (£ n ,r n )} 



(2.4) 



for any monotone sequences {4}neN, {?*n}neN with i < £ n < r n < r and lim ntoo i n = I, lim ntoo r n = 
r. The endpoints of the interval / = (£, r) are absorbing for X(-), that is, we suppose that on {S < oo} 
the process stays after S at the endpoint where it exits. 

We shall suppose that this weak solution has been constructed on some filtered probability space 
(0,5, P), F = {5 r (t)}o<t<oo which satisfies the usual conditions of right-continuity and augmentation 
by P— null sets. Because of uniqueness in the sense of the probability distribution, the state process 
X(-) in this solution has the s trong Markov property . A tho rough exposition and study of such equations 
appears in the monograph by lCherny and Engerbertl (120051) . 



2.1 The Feller test 

A well-known criterion to check whether the diffusion X(-) may explode or not (to wit, whether we hav e 
P(5 < oo) > or F(S < oo) = 0) is Feller's test; see Theorem 5.5.29 in lKaratzas and Shrevel dl99ll) . 
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This criterion relies on the Feller test function defined as 

v(x) := J (J* exp ^-2 jf " f (u)du) ^J dz ) dy , x E I (2.5) 

for some fixed constant c € I. Feller's test then states that P(S = oo) = 1 holds, if and only if 
v(£+) = v(r—) = oo . Moreover, the finiteness or nonfrniteness of v(-) does not depend on the choice 
of the constant c. We shall discuss alternative characterizations below, in Subsection l4.4l 

We observe next that applying Tonelli's theorem, with F(-) = f f(z) dz denoting the anti-derivative 
of the function f(-), yields 



v(x) 



exp(2F(z)) 



j exp (-2F(y)) J dy^j dz. 



Since the function l/s 2 (-) is locally square-integrable thanks to the assumption (12.21 ). and since the anti- 
derivative F(-) is continuous, we have v(l n ) + v{r n ) < oo and, in particular, 

< oo) = 1 (2.6) 



for all n G N; in fact, we even have E[5 n ] < oo, see Proposition 5.5.32 in lKaratzas and Shrevel (|199lh . 
There are several situations in which the Feller test function v(-) can be simplified: 

• If the function s(-) is differentiable and b(-) = as'{-) for some a£R, then 

v(x)= (s 2a -\z) f X S - 2a (y)dy]dz. 



In particular, we have the following two special cases: 
- a = 1/2: 



BW -' '^r^^-sOT^ 1 ; <2 ' 7) 



- a = 1: 



v{x)= j c m dz - 

If b(-) = as(-) for some a 6 R \ {0}, then 

1 /•« 1 - exp(2a(* - g)) 
^ ' 2a J c s 2 (z) 

If b{x) = as{x)/x for all x £ R for some a G R \ {1/2}, then 



" (l) = r^l ,4) dz - <2 - 8) 



In particular, with a = 0, we have 



|^dz. (2.9) 
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From these observations we obtain the following corollary that will be of use later on. 
Proposition 2.1. Conditions for explosions in special cases. Feller's test simplifies in the following 



cases: 



(i) Suppose that the function s(-) is differentiable (without loss of generality, we then assume that 
s(-) > on I), and that b(-) = s'(-)/2 . Then ¥(S = oo) = 1 holds, if and only if for some c £ I 
we have 

i j r i , 

-dz = oo = / —r^dz . 



s(z) ' J c s(z) 

(ii) Suppose that I = (0, oo) and b(-) = 0. Then v(0+) = oo holds, if and only if we have 

jfg^cc; (2.10, 

moreover, v (oo) = oo. Therefore, ¥(S = oo) = 1 holds in this case, if and only if (12.101 ) does. 

(Hi) Suppose that I = (0,oo) and b(x) = s(x)/x for all x £ I. Then v(po) = oo holds, if and only 
if we have 

-dz = oo ; (2-11) 



moreover, v(0+) = oo. Therefore, ¥(S = oo) = 1 holds in this case, if and only if ( 12.111 ) does. 



Proof. Part (i) is an application of the identity in (12.71 ) and Feller's test. For part (ii), we first observe that 
the condition v(0+) = oo , along with the representation in (|2.9t with c = 1, imply that ( 12.101 ) holds. For 
the reverse direction, we assume that (12.101 holds and define m(y) := s~ 2 (z)dz for all y £ [0, 1]. 
If limsup 3/ | (ym(y)) < oo holds, then v(0+) = oo by ( 12.91 ) with c = 1. If limsup^ (?/m(y))) = oo 
holds, we observe that we may rewrite (|2.9t as 

v(0+) = lim / X dz = lim ( / / d y ) dz = m(y)dy > / m(y)dy > t)m(rj) 

xioJ x 5 z (z) xi0\J x J x s 2 (z) ) J J 

for all rj G [0, 1], which yields v(0+) = oo. Moreover, we have v(x) > (x — 2) J^ 2 s _2 (z)dz for all 
x > 2, which yields i>(oo) = oo. 

For part (iii), we use (12.81 ) with a = 1 and c = 1 to obtain the representation 



^^dz = r f i r ^-d^ ) d y . 



In the same manner as above we have v(0+) = oo under (12- lib , and need only show that (12.111 ) implies 
v (oo) = oo . We may assume again limswp y ^ 00 (k(y)/y) = oo, where k(y) = z 2 s~ 2 (z)dz for all 
y £ [1, oo), as otherwise the statement is clear. Under this assumption, we obtain from (12.121 ) that 

f(oo) > [°° \k{z)dz > k(y) \dz = ^ 



holds for all y > 1, which concludes the proof. □ 

Of course, additional statements in the form of the last proposition can be proved; we focused here 
on those we shall use later on. An alternative proof of the condition in ( 12. 101), under slightly stronger as - 
sumptions and using the Ray-Knight theorem, appears in Theorem 1.4 of lDelbaen and Shirakawal (120021) . 
In Example 15.31 we discuss the setup of part (iii) in Proposition 12.11 and its connection to part (ii) will 
then become clearer; see also Corollary [ 
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2.2 A diffusion in natural scale 

Let X°{-) denote the state space in a solution of the equation (12.11 ) without drift, that is, 

x°(-) = S+ I s(X°(t))dW°(t), 
Jo 



(2.13) 



where W°{-) denotes a Brownian motion. Thanks to our assumptions on the function s(-), this equation 
admits a weak solution on some filtered probability space (S7 ,# ,P°), F° = \^ {t)\n<t<no- Th is 



soluti on is unique in the sense of the probability distribution (see Theorem 5.5.7 in lKaratzas and Shreve , 
19911) up until an explosion tim e, and without any absorbing points in the interval / (see Corollary 4.20 in 
Engelbert and Schmido, Il99lh . Of course, the endpoints of the interval are absorbing: once the process 
reaches one of them, it stays there. 

Let S° and 5° denote the sto pping times that are denned as in (12.41 ) with X(-) replaced by X°(-). 
We recall from Theorem 5.5.4 in lKaratzas and Shrevel dl99lh that ¥°(S° = oo) = 1 holds if i = -oo 
and r = oo . 

In order to see the claimed uniqueness in the sense of the probability distribution from first principles, 
and also draw some conclusions that will be valuable later on, we shall also consider a "stopped" version 
of the stochastic integral equation in (12. 13b . namely 



[X(t))dW(t). 



(2.14) 



Here # : C([0, oo)) — y [0, oo] is a given stop-rule denned on the space C([0, oo)) of continuous func- 
tions ro : [0, oo) — > R, with values in the extended real half-line and with the property 



{ro G C([0, oo)) : 0(tt) < t} G B t := ft\B), for all < t < oo. 



(2.15) 



We have denoted by B the Borel a— algebra generated by the o pen sets in C(\0, oo)), and by ipt the 
mapping (ip t tv)(s) := ro(s A t), < s < oo; see Problem 2.4.2 in lKaratzas and Shreve! (|199 lb . We shall 



assume that tv(t) G / holds for all ro G C([0, oo)) and < t < i?(ro). A weak solution of the stochastic 
equation (12.141) consists of a filtered probability space (0,-5", P), F = {5(i)}o<t<oo and of a pair of 
continuous, adapted processes X(-), W(-) on it, such that W(-) is standard Brownian motion and (12.141 ) 
holds. We note that g = is then a stopping time of the filtration F; conversely, we obtain a solution 

of (12.141 ) by stopping the state process X°(-) of any weak solution of the stochastic integral equation 
(12- 13b - Thus, uniqueness in distribution for (12.131 ) will follow once this property has been established for 
the equation (12- 14b - 

Towards this end, let us consider any weak solution of the equation (12.131 ) and denote g = &(X). 
The solvability of (12.141 ) implies that the time change A(-) := f Q s 2 (X(t))l^ s>t jdt is well-defined, 
and we note that this process is the quadratic variation of the continuous local martingale M(-) = 
f s(X(t))l{n^ qdW(t). According to the Dambis-Dubins-Schwarz theory (see Theorem 3.4.6 and 
Problem 3.4.7 in lKaratzas and Shrevel. 1 1 99 ll) there exists a standard Brownian motion B(-) on (an exten- 
sion of) the underlying probability space, such that 



X(.)=e+ / s(X(t))l {e>t} dW(t)=£ + B(A(-)). 
Jo 

We consider now the inverse time change T(0) := mi{t > : A(t) > 9} for < 6 < A(g) and 
r(0) := oo for 9 > A(g), and note that 



r'(9) 



A'(T(9)) s 2 {X(T(9))) ' 



thus 



T(9) 



f 

Jo 



dr 



s 2 (£ + £(r)) 
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for < 9 < A(g). Next, we define the function s*(x) := $(x)lr£ r \(x) + 1 R \ u r \ (x) and the corre- 
sponding time change 



r*(9) 



dr 



5 2 ^ + B(r)) 



< T(9), for all < 9 < oo, 



along with its inverse := inf{<9 > : T*(9) > t} > A(t) for all < t < oo. We have 

A(- A g) = A*(- A g), therefore the process X*(-) := £ + i?(A*(-)) satisfies X (• A g) = X*(- A £>) and 



g = (X(.)J = A g)j = i9 Aj))=* (*,(■)) . 

The continuous, strictly increasing process T*(-) is adapted to the filtration F B = {3 B (9)}g>o generated 
by the Brownian motion B(-), so (A*(t), < t < oo) is a collection of stopping times of ¥ B . This 
implies that the process = £ + -B(A* (•)) and the stopping time g = i9(X !t (-)) are also $ B (oo)- 

measurable. However, we notice that A(-) = A*{- A g), so the processes ^4(-) and X{-) = £ + £?(A(-)) 
are also J B (oo)-measurable. In particular, the distribution of X(-) is detemiined uniquely. 

Remark 2. 1. Construction of a weak solution for (12- 13b - It is clear how to reverse the steps of this 
analysis, and construct a weak solution of (12.131 ): we start with £ G / and a standard Brownian motion 
B(-), and define the first exit time 



r := inf{9 > : f + B(0) g (£, r)}, 



as well as the time change 



T(9) :- 



o 5 2 ( ? + B(r)) 



dr , 



< 9 < t 



with r(#) := oo for all 9 G [r, oo), just as above. Next, we construct the inverse A(-) of T(-), and from 
it the state process X°{-) = X°{- ,£) = £ + B{A{-)). Finally, we construct the Brownian motion 



W°(t) 



0<t<S° = T(r) 



, 5 2 (? + B(r))' 
and check that the pair (W°(-), X°(-)) satisfies (12.131 ) up until the explosion time 

, (r _ )= f d9 f2L(r,y) 



S° = F( 



-dy: 



(2.16) 



here L(r, y) is the local time accumulated up to r by the Brownian motion £+!?(•) at the site y £ M. □ 
2.2.1 A new proof of the Feller test for diffusions in natural scale 



The representation (I2.161 ) for the explosion time 5° and an application of the Ray-Knight theory, as 
for example suggested in Assing and Send (1991), now yield a new pr oof of Feller's condition i n (12.91 ) 
for the diffusion in n a tural scale X°{-)\ in this conn ection, see also Delbaen and Shirakawal (120021) . 
Mijatovic and Urusov (2012), and the related work in iKhoshnevisan et all (|2006l ). More precisely, 
shall see that F°(S° = oo) = 1 holds as follows: 



we 



(i) with 



-oo and r = oo ; 
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(ii) with 



-oo and r < oo , if and only if for all n G N we have 

IT 



y 



dy 



oo ; 



rn s 2 (y) 

(iii) with t > — oo and r = oo , if and only if for all n G N we have 

y 



s 2 (y) 



dy 



oo ; 



(2.17) 



(2.18) 



(iv) with i > -oo and r < oo if and only if both (12.171 ) and (12.181) hold for all n G N . 

The situation of (i) is clear. For the other cases, we recall the local time L(-, •) and stopping time r of 
Remark l2~n Then the function y i— >• L(r, y) is continuous and with compact support, so the local square 
integrability of s (•) gives //; 2L(T,y)s- 2 (y)dy < oo, a.e. Thus, by (T2T61 we shall have 5° = oo 
a.e., if and only if 



2L(T,y) 



and 



2£(r,y) 



dy 



oo 



holds on the event {£ + -B(t 
dy = oo holds on the event |^ + i?(r) 
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•}■ 



(2.19) 



(2.20) 



In t he following, we shall show t he equivalence of (12.171) and (12.201 ) (provided that r < oo) in the 
spirit of lMijatovic and Urusovl (120121) : an analogous treatment then covers the equivalence of (12.181 ) and 
(12.191) (provided that I > — oo). Thus, we now assume, without loss of generality, that r < oo, I = — oo, 
and £ = r n . Under this assumptions, we know from the Ray-Knight theoiy (see Theorem 6.4.7 in 



Karatzas and ShreveU 19911) that the process 2L(r, r — t), < t < r — r n has the same distribution as 
the two-dimensional squared-Bessel process Bf(t) + £? 2 (*)>0 < t < r — r n , where B\(-) and i?2( - ) 
are independent standard Brownian motions. Therefore, with B(-) a standard Brownian motion, (12.201 ) 
is equivalent to the validity of 

- r " B 2 {t) 



21 +\ dt 

5 z (r - t) 

We note that an application of Tonelli's theorem yields 

B 2 {t) 



oo, 



E 



5 2 {r-t) 



dt 



s 2 (r-t) 



a.e. 



dt 



s 2 (y) 



dy. 



Therefore, (12.201 ) implies (I2.17I ). For the reverse implication, assume that (12.201 ) fails; then the event 

■ r » B 2 {t) 



A 



dt < M 



5 2 (r — t 

for some sufficiently large M > 0. Since 



has positive probability 7 := P(A) > , 



E 



B 2 (t)l A > s 2 t ■ p (a n <|V(t) > 5 2 t^ = 5 2 t ■ (p(A) - p (a n {s 2 (i) < 5 2 })) > s 2 t I 



holds for sufficiently large S > 0, Tonelli yields again that 



M > E 



B 2 (t) 
5 2 (r — t 



■dt 



> 



<5 2 7 



5 2 {r-t) 



dt 



2 



r n * 2 (y) 



dy 



holds, namely, that ( 12.171 ) must then fail as well. This shows that ( 12.171 ) implies ( 12.201 ), and concludes 
the argument. 
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3 Transformations of probability measure and space 



To compute probabilistic quantities, for instance distribution functions, one relies on mathematical tools 
such as transformations (distortions) of probability measure, space, or time. As we are interested in the 
distribution of an explosion time, distortions of time do not seem helpful for its computation. In this 
section, we discuss possibly helpful ways to change the probability measure based on the Girsanov-Van 
Schuppen-Wong theorem, and to distort the space based on methodologies developed by Lamperti, Doss 
and Sussmann. 



3.1 A generalized Girsanov theorem 

We di scuss here a generalized version of Girsanov's theorem, which appeared in Section 3.7 of McKear] 
( 1969 1) under condit i ons co nsiderably stronger than those imposed here; see also Exercise 5.5.38 in 
Karatzas and Shrevd (119911) for the special case s(-) = 1 and t = — oo, r = oo. McKean's version can 



be considered a "weak" result, as it provides a distributional identity. iRuil (|2013bl ) uses a related "strong" 
version to provide a proof of the sufficiency of the Novikov and Kazamaki criteria for the martingale 
property of stochastic exponentials. 

First, we recall the finiteness of integral functionals under additional square-integrability assumptions 
on certain related functions. 

Remark 3.1. Finiteness of integral functionals. Let us assume that the mean/variance ratio function f(-) 
is locally square-integrable on /. Furthermore, let us denote by A x (T,y) the local time accumulated 
during the time interval [0, T] by the semimartingale X(-) in (12.11 ) at the site y G /. We may assume that 
the random field A x (- , •) is jointly continuous in bot h its arguments; see for exam ple the statement and 
proof of Theorem 3.7.1 as well as Exercise 3.7 .10 inlKaratzas and Shreve (I199ll) . From the occupation 
time density formula (cf. Theorem 3.7.1 (hi) in lKaratzas and Shrevel Il99ll) and (12.2b . we have 



J b 2 (X(t))dt = ^ (?) {X(t))d(X)(t) = J (fl (in , rn) ) 2 (X(t))d(X)(t) 

= 2 f(y)A x (T,y)dy<2 max (A x (T,y)) • / f(y)dy < oo (3.1) 

Jir, in<V<r n Jo 



-a.e. on the event {S n > T}. 
We are now ready to state and prove a first result. 



□ 



Theorem 3.1. Generalized Girsanov theorem. Suppose that the mean/variance ratio function f(-) is 
locally square-integrable on I. For any given T G (0, oo) and any Borel set A G By(([0, oo))), we have 
then 



(X(-) G A, S > T) = E° 



exp b (X°(t)) dW°(t) -j£* 2 ( x °(t)) dt) 1{X"(.)6A, s°>T} 



In particular, if both diffusions X(-) , X°{-) are non-explosive, i.e., if ¥(S = oo) = P°(S , ° = oo) 
then the exponential P°— local martingale 



(3.2) 
= 1, 



exp ^ b (X°(t)) dW°(t) ~j£b 2 (X°(t)) d*) , < T 



< oo 



is a true P°— martingale. 
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Proof. We fix T G (0, oo) and a Borel set A G j6 t (C([0, oo))) in the notation of (12- 15b - In addition to 
the stopping times of ( 12.41) , we consider the stopping times 



T n := S n A inf jt > : ^ b 2 (X(s))ds > n j , n G N, (3.3) 

as well as stopping times S° and T°, defined in the same manner as in (I2.4t and (13.3b . but with X{-) 
replaced by X°{-). 

We note that 0J} implies {5 > T} = [J neN {T n > T}, modulo P; similarly, we have {S° > 
T} = UneNi^n > modulo P°. In conjunction with the monotone convergence theorem, these 
observations imply that, in order to prove (13.21 , it is sufficient to show 



P(X(-) G A,T n > T) = E° [Z (T)l {X0{ . )6AjTn0>T} ] (3.4) 
for all n G N, where we have set 



:= exp (jf b(X°(t))d^°(t) - ^ jf b 2 (X°(i))dtj 



In the following we shall prove (13.41 ) for fixed n G N. Towards this end, we define the processes 

r-AT n 

W(-) := / b(X(t))dt + W(-), 
Jo 

/ r-AT n i r-AT n \ 

Z(-) := exp Mf b(X(i))dW(i) - - y b 2 (X(t))dfj 

and the P-local martingale 

/ r-ATn i r-ATn \ 

L(-) := exp (- ^ b(X(t))dW(t) - - ^ b 2 (X(i))dtj 



We note that L(-) is a strictly positive martingale; cf. Corollary 3.5.13 in iKaratzas and Shrevel (119911) . 
Therefore, dQ = L( T) dP defines a new proba bility measure Q on (£l,$(T)). Girsanov's theorem 
(see Theorem 3.5.1 in Karatzas and Shreve , ll991 ) yields that W(-) is a Q-Brownian motion. Moreover, 
X(-) := X(- A T n ) is a solution of the stochastic differential equation 

r - AT " — 
*(■)=£+ / *(*(*)) dW(t). 
J o 

This observation and the discussion in Subsection l2.2l implv that the Q-distribution of X(- A T)lr Tn>T \ 
is the same as the P°-distribution of X°(- A T)l{ T o >T }. 
Finally, we note that 

■AT n i r-AT n 



Z(-) = exp (/ f(X(s))dX(s) --] b 2 (X(s))dsj , 

(r-AT° i f-AT" \ 

J " KX°(s))dX°(s) --] " b\X°{s))ds} 



are also nonanticipative functional of X{- A T n ) and A T°), respectively, so we have 

P(X(-) G A, T n > T) = E [L(T) • Z(T)l {x( . )6A|Tn>T} ] 

= E Q [Z(T)l{ X( . )6AiTn>T }] = E° [Z°(T)l{ X o(. )6AjT o >T }] ; 

this yields (13.41 ) and concludes the proof. □ 
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3.1.1 Feynman-Kac representation 

With A = C([0, T]) in (I3T21 ). the distribution of the explosion time 5 in (l2T4l > is given by 



P(5 > T) = E° 



exp b(X°(t))dW°(t) - i jT b 2 (X°(t))di) l {5o>T} 



(3-5) 



if the function f(-) in (12.31) is locally square-integrable on the interval I. Let us assume in this subsection 
that the function f(-) is actually continuous and continuously differentiable on /, and denote by F(-) = 
J f (x)dx its antiderivative. 

On the event {S° > T} we have then 

F(X°(T)) - F(0 = £ f(X°(t))dX°(t) + 1 -£f(X°(t))d(X°)(t) 

= f b(x°(t))dw°(t) + \f T (fV) (x°(t))dt , 

Jo 1 Jo 



so (13.51 ) gives 

P(5 >T) = exp(-F(0) - IE 
with the notation 



exp (f(X°(T)) - V(X (t))dt 



1 {S°>T} 



(3.6) 



V{x) := \ (b 2 (x) + f (x)s 2 (x)) = \ 5 \x) (f 2 (x) + f(x)) = \ (b 2 (x) + b'(x) 5 (x) - b(x) 5 '(x)) 

for all x G I. In other words, the distribution of the explosion time is determined completely by the joint 
distribution of X°(T) and J Q T V(X°(t))dt on the event {5° > T}, for all T G (0, oo). 

Remark 3.2. Non-explosive X°(-). When P°(S , ° = oo) = 1, the expression of (I3.6t takes the simpler 

r-T 



form 



»(5 > T) = exp(-F(0) • E 



exp \F(X°(T)) - J o V(X°(t))dt 



In the special case s(-) = 1 we have X°(-) = £, + W°(-), so finding the distribution of the explosion time 
S as in ( 13.61 ) amounts then to computing the joint distributions of certain Brownian functionals. □ 



3.2 Explosions as Brownian exits via Lamperti-Doss-Sussman-type transformations 

In the previous subsection we saw how to remove drifts by changing the underlying probability measure. 
We discuss now ways to transform the dispersion term into a constant, by distorting the space as Y(-) = 
h(X(-)), for some strictly increasing and continuous function h : (£,r) — > (£,r) and suitable — oo < 

£ < r < oo. 

We shall assume that the function s(-) is continuous and continuously differentiable on the interval 
I = (£, r); without loss of generality, we shall also assume that s(-) is strictly positive. We shall consider 
the function 

h c (x) = [ -Ldz, xe (£,r) (3.7) 

Jc 5 \ z ) 
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for some c G (£,r). We observe that h c (-) is strictly increasing and twice differentiable. We set £ = 
h c (£) := lim^ h c (x) G [—00,00) and r = h c (r) := \\m x ^ r h c (x) G (—00,00] and define the process 
Y(-) via 

Y(t) := h c (X(t)) for all t G [0, 5) and Y(t) = \imh c (X(u)) for all t G [S, 00). 

MIS' 



It is clear that lim^s Y(t) G {£, r } holds on {5 < 00}, and that the new process Y(-) leaves its state 
space / := (£, r) at exactly the time S. In particular, the (distribution of the) explosion time S of X(-) is 
exactly the (distribution of the) explosion time S of Y (•). We recall from part (i) of Proposition 12. II that 
P(5 = 00) = 1 if and only if £ = h c (£) = —00 and h c (r) = r = 00. 

With $ c : I — > I denoting the inverse function of h c , simple stochastic calculus yields that 



dY(t) 



hold for all t G [0, S). In particular, with the function u : I — )■ I defined by 

v( y ):=b(My))-^s'(My)), y^i 

we have the simple dynamics 



dY(t) = u(Y(t))dt + dW(t), Y(0) = h c (0 = e 



for all t G [0, S). As Lamperti ( 1964h stresses (see also Section 3.4 in McKean , 1969 ). this equation 
can be solved pathwise by simple Picard iterations, without any need for stochastic integration or other 
probabilistic tools, as long as the function v(-) is Lipschitz continuous. 

Applying first this transformation of space, and then the change of measure of the previous subsec- 
tion, reduces the problem of computing the distribution of first exit times to a purely Brownian compu- 
tation - although often a hard one, due to the presence of path integrals. Alternatively, a weak additional 
condition on the functions b(-) or v{-), respectively, allows the characterization of the explosion times 
of X(-) as first exit times of a Brownian m o tion from an open (in general , time- dependent) domain via 



the approach pioneered bv Lamperti ( 1964 ). Doss 



stance, in Proposition 5.2.21 in iKaratzas and Shreve 



9771). and ISussmannl (|1978l) and developed, for in 



(119911) . Lamperti, Doss, and Sussmann introduced 



this method to study the pathwise solvability of stochastic differential equations. 

• For completeness, let us sketch the broad outline of this approach. It is checked easily that the map 
(w, c) 1 — y q d c {w) is differentiable on the domain D := {(w, c) : w G Range(/i c ), c G /}, in fact with 





M0)=C, ^ c (0) = l 











w 



d 2 

dw 2 



Mw) = Cas')(0cM)> 



(3.8) 



and 



dc ° 



(w) 



exp 



exp 



*' (My)) ■ &My) 
<My)) 



dy 



«(*eW) S(1? C H) 



<Mo)) 



0(c) • 
(3.9) 
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We assume that for a given continuous function fj, : I — >■ M, there exists almost surely a solution 
C(-) of the integral equation 

C(t) = £ + /% (tf (#(s))) «(C(a))da , < t < S , (3.10) 
■/ o 

where 

S := inf {i > : (W(t),C(t)) $ D} . 

For example, such a solution exists if the function /x(-) is locally Lipschitz continuous on /. If we define 
now 

X(t) :=0 O (t)(W(*)) for all < t < 5, 
it is checked easily that X(-) solves (127Tb with b(-) = s'{-)/2 + jLt(-); to wit, 

dX(i) = s(X(t)) (dW(t) + (ys' + ^lWjdtj , X(0)=£gJ. (3.11) 

This follows from elementary stochastic calculus, using (T3.8I) - (13.10b . 

The approach just outlined resembles mostly that of iDossI (119771) : however, our representation of 
the drift lets the ordinary differential equation in (B.lOt look slightly easier than in that paper. We also 
observe that X(-) does not depend on the path W(-) through the subindex C(-), if and only if fi(-) is a 
real constant. 



4 Analytic properties of the explosion time distribution 

In this section, we shall discuss analytic properties of the function U : (0, oo) x / — > [0, 1], defined via 

U(T,0 :=P 5 (5 > T) (4.1) 

for all (T, £) E (0, oo) x /. Here and in what follows, we index the probability measure by the common 
starting position £ € I of the diffusions X(-) and X°(-). 



4.1 Continuity 

The question of continuity of the function U(- ,■) is of interest in itself; it also will be important for our 
arguments later on. Since 1 — U(-, £) is a distribution function, it is right-continuous for all £ G /. In this 
subsection we shall see, without any further assumptions on the coefficients s(-) and b(-) beyond those 
imposed in Section 1, that {/(-,•) is actually jointly continuous in its two arguments. 

We start with a technical result in Lemma 14.11 In particular, the property in (I4.4t is well known for 
re gular, one-dimens i onal d iffusions; it is discussed, for instance, in the "matching numbers" Section 3.3 
of lltoand McKearJ Jl965h . It is not hard to prove from first principles, so we present here a simple 



argument. We then prove in Lemma l4~2l the continuity of the function U(- , ■) in the first component, as 
a function of time only. Finally, in Proposition 14 . 3 1 we establish the joint continuity of £/(•, •) based on a 
simple "coupling" argument. 

Lemma 4.1. Diffusions hit nearby points fast. With the stopping times 

H x := inf {t G [0, oo) : X(t) = x} , x G I, (4.2) 
H XiV := inf {t G [H x ,oo) : X(t) = y} , (x,y) G I 2 , 
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for any given e > there exist x\ = x\(e) G (— oo, £) and X2 = X2(e) G (£, 00) such that 



P^H x ^<e) >l-e, z = 1,2. (4.3) 
In particular, for all 5 > 0, we have 

lim Py (He <S) = 1. (4.4) 
Proof. We first show that we have P^-A) = for the event 



i:= weSJ:] i?M > such that min X(t, u) > O ■ (4-5) 

1 te[0,ii(u/)] J 

Without loss of generality, we may assume that < min^o^^)] Z(t, lo) < max^o^^)] Z(t,ui) < 2 
for all uj G $7, where Z(-) is the stochastic exponential of f Q b(X(t))dW(t). Thus, by an equivalent 
change of measure argument, we also may assume that b(-) = 0. The path properties of standard 
Brownian motion, in conjunction with the observations of Remark l2~T1 a nd the fact that A(t) > fo r all 



t > 0, as defined there, then let us conclude; see also Lemma 2.9.4(h) in lKaratzas and Shrevd (119911) . 

The continuity (from below) of the probability measure P^ then yields the existence of y\ G (—00, £) 
such that Pe(H yi < e) > 1 — e/2. Replacing the minimum by a maximum in ( 14.51 ) and repeating the 
argument, we obtain the existence of X2 G (£, 00) such that P^(H X2 < H yi < e) > 1 — e holds; this 
then implies ( 14.31 ) for i = 2. The existence of the claimed x\ G (—00, £) is argued in the same manner. 
Finally, the strong Markov property of the diffusion X(-) implies 

H, < e) > P f (h, f < e 



r y <v t) ^ irt yny£ 

that is, the probability of the event that the diffusion X(-) started at y hits £ before time e dominates 
the probability of the event that X(-) completes a round-trip from ^ to y and then back to £, before time 
e. We now fix 5 > 0,e G (0,5) and the corresponding x\ G (— 00, £) and X2 G (£,00). Then for all 
y G (xi, X2), applying (I4.3I ). we have 

P y (H{ <S)>F y (H^<e)>l- e, 

which proves (14.4b - □ 

Lemma 4.2. Continuity of U(-,-) as a function of time. The function T \-t U(T,£) is continuous on 
[0, 00), for any given £ G I. 

Proof. We fix (T, £) G (0, 00) x / and observe that it is sufficient to show p := P^(S = T) = , due to 
the right-continuity of the function U(- , £). Let us consider any strictly increasing sequence of stopping 
times = 

jy(o) < < #(2) < . . . _ We then have? 

again by the strong Markov property of the 

diffusion X(-), 

1 > P 5 (5 G {r + ff« : i G No}) = ]T P ? (s - # (i) = r) = ^ E $ [p x( „ W) (5 = T) 

(x(ff«) = e). 

Thus, in order to show the statement, it is sufficient to construct a strictly increasing sequence of stopping 
times sucn that P^ (X(H^) = £) does not converge to zero as % increases. We shall construct 
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such a sequence inductively, by "stitching together" the round trips of Lemma 14.11 Towards this end, 
consider a sequence {(?j}ieN C I such that (14.31 ) holds with x\ replaced by qi and e replaced by 1/2*. 
Next, we define the stopping times 

:=inf{iG [tf^oo) : X(t) = ft} , 
#W := inf {t € [H^ ,oo) : X(t) = . 

Then using the Markov property of the diffusion X(-), and conditioning on the event = £}, 

we obtain 

n (|>K) 



> ... > 



which does not converge to zero as i increases. This concludes the proof. □ 

Alternatively, we could have proved the previous lemma via the methods developed in Section [3] 
under the additional assumption of continuous differentiability of the function s(-), for instance. We 
would first simplify the stochastic differential equation to one with s(-) = 1 via the Lamperti-Doss- 
Sussman-type transformation. In a second step, we would find a suitable localization and remove the 
drift by a change of the probability measure. We then could conclude by using properties of Brownian 
motion. The advantage of the proof we provided, is that it does not require any assumptions on the 
function s(-). 

Proposition 4.3. Joint continuity of £/(•,•). The function (T, £) h-» U(T,£) is jointly continuous on 
[0, oo) x I. 

Proof. We fix some e > 0, a pair (T, £) G [0, oo) x /, and a sequence {(t n , £ n )}neN C [0, oo) x / such 
matlim ntoo (t n ,£ n ) = (T,f). 

We start with the case T = 0. We need to show that lim ntoo F^ n (S > t n ) = 1. With 5 = 
S(£) A S(r), the minimum of the explosion times of X(-) to £ and r, respectively, let us choose some 
7] G (€, £) and observe that, for sufficiently large n G N, we have the upper bound 

P*„0SW > t n ) > F V (S(£) > t n ) = U(t n , v ) 

which converges to U(0,?]) = 1 as n f oo due to the right-continuity of the function U(-,rj). Thus, 
\im. n ^ O0 'F^ n {S{l) > t n ) = 1; similarly lim^oo Pg n (5(r) > t n ) = 1, and this proves the claim for 
T = 0. 

We now assume T > 0. From Lemmagl there exists 5 G (0, T/2) so that \U(t, £) - E/(T, £)| < e 
holds for all t G (T - 25, T + 25). Without loss of generality, we assume \t n - T\ < 5 for all n G N. 
Next, we observe that the strong Markov property of X(-) implies 

\U(t n ,Cn) - U(t,£)\ = |e s „ [i {s>tn} - u(t,o] I 



< E* 



\u{t n -h^o- U(t,0\l{H 6 < S } +p 5 „ > s) 



<e + F^(H^>5); 

here i/g , defined as in (14.2b . is the first hitting time of £ by the process X(-). Letting n tend to infinity 
and applying (14.41 ) concludes the proof. □ 
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4.2 Connections with parabolic partial differential equations 



In this section, we study conditions implying that the function U(-,-) solves the Cauchy problem for the 
linear, parabolic partial differential equation 

5 2 (x) d 2 u 



du 



dlA 

dx2 -(r,x) + b(x)s(x) — ( T ,x) 



with an appropriate initial condition, namely 

U(0,x) 

We start with an existence result. 



x G I. 



(4.6) 



(4.7) 



Lemma 4.4. Existence of a classical solution. Assume that the functions s(-) and b(-) are locally 
uniformly Holder-continuous on I. Then for any bounded, continuous function g : (0, oo) x I — > oo and 
any n G N, the Cauchy problem of (14.61 ) has a unique classical solution IA{- , •) of class C 1,2 ((1 /n, oo) x 
(Ai> r n))» subject to the boundary conditions 

U(l/n,x) = g(l/n,x), U{r,i n ) = g{r,t n ), and U(r,r n ) = g{r,r n ) 

for all (t,x) G (1/n, oo) x (l n ,r n ). Moreoever, the solution U{-,-) is bounded on [l/ra,oo) x [£ n ,r n ]. 



Proof. The continuity of the function s(-) yields mm x& y n ^ \s(x)\ > 0. Moreover, our assumptions 
imply that the functions s 2 (-) and b(-)s(-) are als o uniformly Hold er-continuous on [^n,r n ]. Thus, the 
existence and uniqueness result of Theorem 3.9 in iFriedmanl (| 1 9641) and the maximum principle of The- 
orem 2. 1 in this same book yield the statement. □ 

We can now show that the function U (•, •) of d4.il ) solves the Cauchy problem of ( 14.61 ), (I4.7I ). 

Proposition 4.5. Stochastic representation of the solution to the Cauchy problem of (14.61 ). (14.71 ). 

Under the assumptions of Lemma \4.4\ the function U(- , •) is of class C([0, oo) x I) n C 1,2 ((0, oo) x /) 
and solves the Cauchy problem of (14.61) . (14.7I ). 

Proof. We have shown the continuity of the function U(-, •) in Proposition 14.31 We now fix (T, £) G 
(1/n, oo) x (£ n ,r n ) for some n G N and show that the function U(- , •) satisfies the Cauchy problem of 
(14.61 ) in (1/n, oo) x (£ n ,r n ), which then yields the statement. 

Applying Lemma 1441 with g(-, •) = U (•, •), we see that this Cauchy problem has a bounded classical 
solution U{-, ■). Simple stochastic calculus then implies that U(T — (t A p), X(t A p)), < t < T is a 
bounded P^-martingale, where p denotes the smaller of T — 1/n and S n . Optional sampling gives 

U{T,Q=Ez[U(T-p,X(p))] = E i [U(T-p,X(p))] =U(T,£), 

where the last equality is a consequence of the strong Markov property of the diffusion X(-). We have 
shown that U(-, •) coincides with U(-, •), and thus solves the Cauchy problem of (I4.6l l. (I4.7I ). □ 



The proof of Propositio n |4^5Jresembles the argu ments inUanson and Tvskl (I2006I). For sim i lar res ults, 



Heat h and Schweiz er (200Q) and 



see Section 3.5 in lMcKeanI jl969l ). Theorem 5.6.1 in lFriedmanI dl976l) . 

Ruld2013al \ We emphasize that the Holder coefficient in Lemma l4!4l need not be 1/2, as often postulated 
in related questions. 

This Cauchy problem of (14.6b . (14.71 ) admits the trivial solution U{-, •) = 1; it may have lots of 
other solutions. The one we are interested in, the function U(-,-) defined in (14.1b . turns out to be its 
minimal nonnegative sol ution. The follo wing characterization of this function is analogous to the re- 



Fernholz and Karatzasl (1201 (J . 1201 ll) and 



suits in Problem 3.5.1 of McKe arj 1 1969) and in Exe rcise 4.4.7 of iKaratzas and Shrevel (|l99ll) : see also 



Ruj (1201 3ah . 
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Proposition 4.6. Upper bounds on {/(-, •), and minimality. The function U {•,•), defined in (14.11 ) as 

the tail of the distribution function of the explosion time S, is dominated by every nonnegative classical 
supersolution U(-,-) of the Cauchy problem of (14.61) . (14.7b . 

Whereas, under the conditions o f Lemma \4.4\ the function U(- , •) is the smallest nonnegative classi- 
cal (super)solution of (14.6b . (14.7b - 



Proof. Consider any continuous function U : [0, oo) x I — > [0, oo) of class C 1,2 ((0, oo) x I) that satisfies 
the partial differential inequality 

du . , 5 2 { x )d 2 u r . t . w .aw. , , . , . r 

^t(t,x) > — ^— (t,x) + b(x)5(x) — (r,x); (t,x) g (0,oo) x J, 

as well as the initial condition 14(0, •) > 1. For any given T E (0, oo), it is checked readily on the 
strength of this inequality that the process U(T — (t A S n ),X(t A S n )) for < t < T is a local 
Pg— supermartingale; as it is nonnegative, this process is actually a true P^— supermartingale, so 

U(T,Z)>Es[U(T-(TAS n ),X(TAS n ))] > E ? [l {Sn>T} U(0,X(T))] > P^ (S n > T) (4.8) 

by optional sampling. Letting n f oo, we obtain the first claimed result from monotone convergence; the 
second then follows from Proposition 14.51 □ 



If there are no explosions, then the Cauchy problem of (I4.6I ). (14.7b actually has a unique bounded 
classical solution. 

Proposition 4.7. Unique bounded solution. Assume that the function U(-, ■) = 1 is the smallest non- 
negative classical solution of the Cauchy problem (14.6b . (14.7b . and let V(-, •) be any bounded classical 
solution of this Cauchy problem. Then we have V(-, •) = £/(•, •) = 1. 

Proof. Let V : [0, oo) x I — > [-K, K] be a bounded classical solution of the Cauchy problem (14.6b . (14.7b . 
for some K £ (0, oo). Then the function V(-, •) = (V(-, •) + K) /(l + K) is a bounded, nonnegative 
classical solution of the same Cauchy problem, so by assumption we must have V(-, ■) > 1. 

Let us assume that V(- , •) is not identically equal to the constant function U(-,-) = 1, so that we 
haveifa := sup (irJ , )g[0oo)x(Ar) V(t,x) G (1,2). The function V"(- , •) = (K 2 -V (;•))/ (K 2 - 1) is 
then a classical solution of the Cauchy problem (14.6b . (14.7b with values in [0, 1], and not identically equal 
to the constant function U(- , ■) = 1; but this contradicts the assumption that the function U(- , •) = 1 is 
the smallest nonnegative classical solution of the Cauchy problem (14.6b . (14.7b . □ 



We note that it is not po ssible to remove the boundedness assumption in Proposition I4.7[ see for 



example iRosenblooml (1 19580 . 



4.3 Connections with second-order ordinary differential equations 

Let us consider now, for any given real number A > 0, the Laplace transform or "resolvent" of the 
function U(-, £) in (14.1b . namely 



roc poo i . 

/ exp(-\T)U(T,£)dT = exp(— AT)P^(5 > T)dT = — (1 — Eg [exp(— AS*)] 
Jo Jo * ^ 



(4.9) 
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Proposition 4.8. Stochastic representation of the solution to an ordinary differential equation. If 

the functions b(-) and s(-) are continuous on I, then the function U\(-) is of class C 2 (I) and satisfies 
the second-order ordinary differential equation 



u"(x) + b(x)s(x)u'{x) - Xu(x) + 1 = 0, x G I. 



(4.10) 



Proof. For some fixed n G N, we consider the ordinary differential equation 

5 2 (x) 

— j-±v (x) + b(x)s(x)v (x) - Xv(x) = 0, x G (£ n ,r n ) 



(4.11) 



with boundary condition v(£ n ) = v(r n ) = 0, and note that it has the unique solution u(-) = 0. To 
see why, let us assume that ( 14.1 11 ) has a non-constant solution v(-), and try to arrive at a contradiction. 
This solution v(-) must have a local maximum or minimum at some y G (£ n ,r n ) with v'{y) = 0; 
assuming that y is the location of a positive local maximum with v(y) > 0, we obtain the absurdity 
> s 2 (y) v"(y) = 2X v(y). This yields the asserted uniqueness. 

Thereom 12.3.1 in iHartmanl (11982b shows now that the differential equation in ( 14.101 ) has a unique 
solution in (£ n , r n ) with boundary conditions u(£ n ) = g\ and u(r n ) = g2 for all n G N and g±, gi G M. 

Finally, we fix a £ G / and a sufficiently large n G N so that £ G (£ n , r n ) and let u(-) denote the 
solution of the differential equation in (14.101 ) with boundary conditions u(£ n ) = U\(£ n ) and u(r n ) = 
U\(r n ). Simple stochastic calculus shows that the process 



M(t) := exp(-A(t A S n )) ( u{X(t A S n )) - j 



< t < oo 



(4.12) 



is a Pg-local martingale; we conclude that M(-) is a uniformly integrable martingale, as it is bounded. 
We recall P^(S n < oo) = 1 from (1231 ). and obtain that 



u(0 



1 

T 



E f 



exp(-XS n ) u(X(S n )) 



1 



E 5 [exp(-A5 n ) • E x{Sn) [exp(-A5)]] 



exp(— A5- 
1 



U x (X(S n )) 



x E 5 [ex P (-A5)] = U X (0 - - 



the result now follows. 



□ 



For related results, see Theorem 5.9.3 inllt6l (]2006|) . and Theorem 13.16 on page 51 in Volume II of 



Dvnkinl (119651) . 



Once again, the ordinary differential equation in (I4.10t may have lots of classical solutions, in addi- 
tion to the obvious u(-) = 1/A. The function of (14.91 ) we are interested in, turns out to be its smallest 
nonnegative classical supersolution. 

Proposition 4.9. Upper bounds on U\(-), and minimality. The function U\(-), defined in ( 14.91 ) as the 

Laplace transform of the tail of the distribution function of the explosion time S, is dominated by every 
nonnegative classical supersolution of the second-order equation (14.1 01 ). 

Whereas, under the conditions of Proposition l?T8l the function U\{-) is the smallest nonnegative 
classical (super)solution of (14.1 01 ). 

Proof. Consider any function u : I — > [0, oo) of class C 2 (I) that satisfies 



5 2 (x) 

— jpV'(s) + b(x)s(x)u' (x) - Xu(x) + 1 < 0, 



x el. 
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Simple stochastic calculus shows that the process M(-), defined in (14- 12b . is now a supermartingale, so 



u(£) > E c 



exp(-AS n ) u(X(S n )) 



/•oo 

/ exp(-At)l {Sn>t} dt 
Jo 



exp(-At)P f (5„ > *)dt. 



Letting n tend to infinity, we conclude that u(£) > E^fJ^ 00 exp(— Xt)¥(S > t)dt] = U\(£) holds for any 
given initial position £ G I; the first claim follows. The second is a consequence of Proposition 14.81 □ 



By analogy with the situation in Subsection l4.2[ in the absence of explosions the second-order equa- 
tion ( 14.101 ) actually has a unique bounded solution. 

Proposition 4.10. Unique bounded solution. Assume that the function U\{-) = 1/A is the minimal 
nonnegative classical solution of the second-order equation (14.101 ). and let u(-) be any bounded classical 
solution of (14. 10b . Then we have u(-) = U\(-). 



Proof. This result can be proved in exactly the same manner as Proposition [477] 



□ 



We note again that it is not possible to remove the boundedness assumption in Proposition 14.101 For 
example, with I = R, s(-) = \/2, b(-) = 0, and A = 1 in (14.101) . the (smallest nonnegative and) unique 
bounded solution is of course U\ (•) = 1, but there is a host of unbounded solutions u K (x) = l+«exp(x) 
for all k G R (nonnegative, as long as k > 0). 

Finally, we remark that solving the ordinary differential equation in (14.1 II) has been the standard way 
to compute Laplace transformations of the distributions of hitting times, a com putation that is a special 



case o f the computa tio n of the distributions of explo sion times; see for example iBarndorff -Nielsen et al. 
dl978b . lKentl <ll978h . or lSalminen and Valloisl h(X)9h . 



4.4 Equivalent formulations for the Feller test 

The next theorem summarizes several of our previous observations. 

Theorem 4.11. Characterization of explosions. The following conditions are equivalent: 

(i) the diffusion process X(-) of (12.11 ) has no explosions, i.e., P(5 = oo) = 1 ; 

(ii) v(£+) = v(r—) = oo hold for the "Feller test" function defined in (12.51) for some c G (£, r). 
If the function f(-) is locally square-integrable on I, then conditions (i)-(ii) are equivalent to: 

( Hi) the truncated exponential ^-supermartingale 

fT 1 r-T 



Z b (T) = exp^ b(X°(t))dW°(t) - ijf b\X° 



(t))dt 



l {S°>T}i 



< T < oo, 



appearing in (I3.5I ). is a P° -martingale. 
If the functions s(-) and b(-) are continuous on I, conditions (i)-(iii) are equivalent to: 



(iv) the smallest nonnegative classical solution of the second-order differential equation in (I4.10I ) is 
u(-) = l/A; 



(iv)' the unique bounded classical solution of the equation in (I4.10I ) is u(-) = 1/A. 
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Whereas if, in addition, the functions s(-) and b(-) are locally uniformly Holder-continuous on I, condi- 
tions (i)—(iv) are equivalent to: 



(v) the smallest nonnegative classical solution of the Cauchy problem (14.6I ). (14.71 ) is IA(-,-) = 1 ; 

(v)' the unique bounded classical solution of the Cauchy problem (14.61) . (14.71) is £/(•,•) = 1. 

Proof. The equival ence of (i) and (ii) is the su bject of the Feller test for explosions (see for example 
Theorem 5.5.29 in iKaratzas and Shrevei Il99ll) . The equivalence of (i) and (iii) follows from Theo- 
rem [3j] Under the stated conditions, the equivalence of (i) and (iv) is covered by Propositions 14. 8 1 and 
|4.9t whereas the equivalence of (i) and (v) is covered by Propositions 14.51 and 14.61 The equivalence of 
(iv) and (iv)' (resp., (v) and (v)') is the subject of Proposition 14 . 1 01 (resp . . Proposition 14.7b - □ 

It would be of some interest to have a more "circular" proof of this result, in particular, a derivation 
of the Feller test (ii) from one of the minimality properties (iv) and (v). 



5 Examples 

Let us consider some examples. 

Example 5.1. Reciprocal of Brownian motion. Let us take / = (0, oo) and 

s(x) = —x 2 , b(x) = —x, 



(5.1) 



implying f(x) = l/x. With a given initial condition £ £ (0, oo), the driftless equation of (12.131 ) takes 
the form 



{X°(t)) 2 dW°(t), 



and is easily seen to take values in / = (0, oo) for all times, as it is identified with the reciprocal 
X°(-) = 1/R{-) of the three-dimensional Bessel process 



dR(t) 



1 



R(t) 



dt + dW°(t), 



R(0) 



(5.2) 



In particular , the condition in (12.101 ) is satisfied, and we have P°(S° = oo) = 1. It has been known since 



the work of 



Johnson and Helms! (11963b that X°(-) is a strict local martingal e; and the connection wit h 



the Bessel process allows the computation of the distribution of X°(T) as in Ekstrom and Tysk d2009h . 
namely 



P(X°(T) e dy) 



exp 



«l/y) ~ (VP) 2 
2T 



exp 



((1/y) + (l/Qf 
2T 



On the other hand, with the choices of (15.ll ). the general equation of (12.11 ) becomes 



X 2 (t)dW(t)+ / X 3 (t)dt, 



(5.3) 



whereas the functions of S ubsection 13 . 1 . 1 1 take the form F(x) = log(x) and V(-) = 0. The condition 
(12.1 II) clearly fails in this case, so we have P(S = oo) < 1 ; in fact, it follows from (13.61 ) that 



P(5 > T) = — ■ E° 



R(T) 



exp 



dr=:U(T,£), < T < oo , (5.4) 
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so in fact P(S = oo) = 0. This is also quite straightforward to check directly, as follows: we observe 
that (15.31 ) is of the special form 

dX(t)=s(X(t))(dW(t) + ^s , (X(t))dt), X(0)=£, (5.5) 

or equivalently (15.171 ) with /j(-) = 0. Following the procedure of Subsection 13.21 we see that (15.31) can 
be solved "pathwise" in terms of the function <&^(w) = (w + (1/£)) _1 , namely as 

Thus, the explosion time S is the first hitting time of the level (— 1/£) by a standard Brownian motion 
started at the origin, that is, the right-hand side of (I5.4I ). Note that the process explodes at S = lim^oo S' n 
where S' n = ini{t > : X(t) > n}; the left endpoint t = of the state space is inaccessible by X(-). 
This is consistent with the observation made in part (i) of Proposition 12.11 since |s _1 (2:)|dz = oo and 
f™\s- 1 (z)\dz = 1 < oo. 

Here, it is easy to verify "by hand" that the function U(-, •), defined in (15.41) . satisfies both the linear 
parabolic equation of (14.61 ). now in the form 



dU . . x 4 d 2 U . . o dU . 

t^t( t > x ) = ( T ' x ) + x -^:( T > x )i ( T > x ) G (°i °°) x ^ 

as well as the initial condition U(0+,x) = 1 for all x £ I. (It also satisfies the lateral condition 
U(T, 0+) = 1 for all T G (0, oo) but this is immaterial, as the left endpoint t = of the state space 
is inaccessible.) The function of (15.41) is not the only classical solution of this Cauchy problem, as 
W(v) = 1 is clearly a solution; from Proposition 14.61 however, U(-,-) is its smallest nonnegative 
classical solution. 

• Let us consider next the second-order ordinary differential equation in (14.101) . written here as 

4 

—u"{x) + x 3 u'{x) - Xu(x) + 1 = 0, x & I. (5.6) 

It is easy to see that a general solution of this differential equation takes the form 

t \ a f 2X \ v f 2X \ 1 
u[x) = A exp I I + B exp I — I + — , xGi 

for some real constants A, B. We are interested in the smallest nonnegative solution U\{-) in (15.6I ). which 
we obtain by first setting B = and then A = — 1/A. We thus have 

&C-)-i(i-«p(-£) 

which is clearly smaller than the constant function 1/A. This illustrates the validity of the characterization 
of an explosive diffusion in (v) of Theorem 14.1 II □ 

Example 5.2. A generalization of Example 15. 11 On the state space / = (0, oo) and with a constant 
k G [1/2, oo), let us consider the stochastic equation 

X{-)=i-f (X(t)) 2 dW(t) + K [ {X(t)) 3 dt. 
Jo Jo 
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This corresponds to s(x) = —x 2 and b(x) = —kx, thus 



f(x) 



F(x) = log(x K ), V(x) 



ax 



with a := k(k — 1), in the notation of Subsection 13.1.11 With the three-dimensional Bessel process 
R(-) = 1/ X°{-) of ( 15.21 ). the representation (13.61 ) now provides the distribution of the explosion time S 
ofX(-) with v := k- 1/2 as 



P(5 > T) = E° 



E 6 



1 



(£J2(T)) 
1 



T dt 

R(t) 2 

v-l/2 



— .(^(T))-^exp 



^-1/4 



1 



(ei?(r)) 



2k 



-1 



l-u 

x exp 



-x 



2T 



w dx - 



in terms of the modified Bessel function of the second type 

{u/2Y +2n 



I„(u) := ^ 



n6N 



n\T{n + 



(5.7) 



(5.8) 



and the Gamma function T(-). Here, we used Q" to denote the probability measure under which R(-) 
is a Bessel pro cess in dimension 5 = 2{v + 1) = 2k + 1 > 2. For this change of measure, see Exer- 
cise XL 1.22 in Revuz and Yoj dl999); the dens ity function p u (-) of the Bessel process R(-) under is 
provided on page 446 of iRevuz and Yon ( 19991) . Alternatively, we could have used Formula (1.20.7) on 
page 448 of Borodin and Salminenl (|2002l ) for the case k > 1. We shall further simplify these computa- 
tions in Example l5.8l 

From Propositions 14 . 5 1 and 14. 6 1 the function in (I5.7t is the smallest nonnegative solution of the linear 
parabolic equation 



dU 



x 



d 2 u 



,dU 



IT,X 



2 dx 2 ^ x) + KX "^ {T ^ 



(t,x) £ (0,oo) x I, 



subject to the initial condition U(0+, x) = 1 for x E I. 



□ 



Example 5.3. Another generalization of Example 15.11 and /i-transforms. The special case / = 
(0, oo) and f(x) = 1/x, or equivalently b(x) = s(x)/x, corresponds exactly to the situation in which 
the dif fusion X(-) is th e /i-transform of the nonne gative local martingale X°(-); see for example iDoob 
or Section 3.3 in lPerkowski and Run (12012b . We observe that Example |5.1| is a special case of this 
setup. The above choices lead to F(x) = log(x) and V(-) = in the notation of Subsection 13. 1.11 and 
(13.61) becomes 



P(5 > T) = j E° [X°(T) ■ l {S o >T }] =j^° [X°{T A 5°)] 



< T < oo 



(5-9) 



for the diffusion X°(-) in natural scale of (12- 13b - This can be seen from first principles, as X°(-) is a 
nonnegative local martingale, which gets absorbed at the origin the first time it reaches it. Consequently, 
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we have F(S = oo) = 1 if and only if the stopped diffusion in natural scale X°(-) of (12.131 ) is a 
martingale. 

Let us recall from part (iii) in Proposition 12. 1 1 that we have v(oo) = oo and P(5 = oo) = 1 , if and 
only if the con dition in (12.11b holds. That is, X°(-) is a martingale, if and only if the condition in (12.1 II ) 
holds; see also lDelbaen and Shirakawal (120021) . Alternatively, we observe that the process Y(-) = \/X(-) 
satisfies the equation 

Y{ ^j-L Y ^-im) dW(t) 

and reaches the origin if and only if 



dz = oo . 



by (12- 1 Ob - This, however, is again equivalent to (12 -lib - 

The equation for Y(-) can be solved by the familial - method of time-changing a standard Brownian 
motion B(-), namely 

i r l 

Y(t) = — + B(A(t)); here A(t) = h 2 (Y(s))ds, < t < oo , 
? Jo 

with h(x) := x 2 5(l/x), is the inverse of the continuous, strictly increasing, real-valued function 

T(u) := J h~ 2 Q + B(vU dv, < u < oo 

with r(oo) = oo. The explosion time S of X(-) is thus related to the first hitting time 

t := inf ju > : B(u) = -j\ 

for the auxiliary Brownian motion B(-) via 

P(S >T)= P(t > A(T)) = P(T(t) > T). (5.10) 

We now consider the special case s(x) = kx p for all x > for some real numbers k > , p > . 
First, by (12.10b . we observe that X°{-) is a martingale and P(S = oo) = 1, if and only if p < 1, that is, 
the function s(-) grows at most linearly. We note that in the case p = 1/2 we have 

dX(t) = n 2 dt + Ky /X(t) dW(t), X(0)=£. (5.11) 

In particular, with k = 2 we see that X(-) is the square of a Bessel process in dimension 5 = 4. 
For the case p > 1 = k it would be very nice to compute the distribution of the random variable 

F{t) = lo il + B(V) ) P ^ 

and thus the distribution of the explosion time 5 via (15.10b . but we do not yet know how to do this in 
general, except in two special cases. 

The first exception is the value p = 2. With this choice h(-) = 1, T(u) = u and P(5 > T) = P(r > 
T) is given by the expression of d5.4t . 
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The second exception is the value p = 3/2, for which iBorodin and Salminenl (120021) provide in 
Formula (2. 19.2) on page 208 the distribution of the random variable 

or equivalently the distribution of the explosion time 

F{S >T) =F{T{t) >T) = l-expf-^J = j ■ E° [X°{T)} , < T < oo (5.12) 



for the diffusion 



*(.)=£+ / (J(i)) 3/2 dW(f)+/ (X(t)) 2 dt ; 
io Jo 



thus = £ + r o (X°(t)) 3 / 2 dW°(t) . We shall discuss a related diffusion in Example[ 

It is now checked "by hand", that the function U(T, £) = 1 - exp(-2/(£T)) in (15.121 ) satisfies the 
linear parabolic equation 

§(T,0 = ^ {T ,0 + e^(T,e , (T,0 € (0,oo) x /, 

subject to the initial condition U(0+, •) = 1 on /; from Propositions 14.51 and 1431 it is the smallest 
nonnegative (super)solution of this equation. □ 



In t his last example, we proved the following result, which appeared under slightly stronger assump- 



tions m 



Delbaen and Shirakawal (120021) . 



Corollary 5.4. Martingale property of nonnegative diffusions in natural scale. Suppose that the 
function l/s(-) is locally square-integrable on I = (0, oo), and that X°(-) is a nonnegative ¥°-local 
martingale which satisfies the stochastic differential equation 



X°(-) = £+ / s(X°(t))dW°(t) 
Jo 



and becomes absorbed at the origin the first time it gets there. 

Under these conditions, X°(-) is a true F c '-martingale if and only if the condition in (12.1 II ) holds. 



This statement has been shown in Example [53J see for example ( I5.9I ). 

Example 5.5. Lamperti-Doss-Sussmann for affinc variance diffusions. Let us consider a modification 
of the stochastic differential equation (15- lib , namely / = (0, oo) and 

K 2 



dX(t) = —dt + Ky/X[t) dW(t), X(0) = £ € (0,oo) (5.13) 

for some k > 0. This corresponds to s{x) = K^fx, b(x) = k/(4-*/x), thus 

1 1 3k 2 

f(x) = — , F(x) = -\os(x), and V(x) = . 

n ' Ax K ' 4 &v ' K ' 32x 
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We note that b(x) = s'(x)/2 holds, so once again the equation (15.131 ) is of the Lamperti-Doss- 
Sussmann form (15.51 ) and can thus be solved pathwise by following the procedure outlined in Subsec- 
tion 13.21 The result is X(-) = (W(-)), now in terms of the function 



< t < S. 



namely X{t) = ( — W(t) + V£ 



It follows that S is the first hitting time of the point (— 2y/£/n) by standard Brownian motion started at 
the origin. Thus, the process X(-) explodes by hitting the left endpoint of the state space; the right 
endpoint oo is inaccessible by the diffusion X(-); that is, X(S-) = 0. Of course, this is again consistent 
with part (i) of Proposition 12. 1 1 since J s~ 1 (z)dz = 2k < oo and s~ 1 (z)dz = oo. 
From (13.6I ). we conclude that 



P(5 > T) = 2 [ 
Jo 



X°(T) 



1 



exp 



2vr 

1/4 

exp 



3^2 

"32" 



T 



(5.14) 



dt 



X°(t) 



1 {S°>T} 



here X°{-) is a martingale and diffusion in natural scale with 

dX°{t) = K^fXHf)dW {t), X°(0) = £. 

Once again, it is easy to check by direct computation that the function of (15.141 ) solves the linear 
parabolic equation of (I4.6I ). namely 



dU 



d 2 u, 



- 2 du 



IT,X 



subject to the initial condition U(0+, x) = 1 for x E I (and to the lateral condition U(T, 0+) = for 
T G (0, 00)). Arguing then by analogy with Proposition 14.61 it is checked that the function of (15.141 ) is 
the smallest nonnegative classical solution of this initial / boundary value problem. □ 

Example 5.6. Diffusion with quartic variance function. Let us consider, for some given real number 
v G M, the stochastic differential equation 

dX(t) = (l + X 2 (t))(dW(t) + (v + X(t))dt) , X(0)=£ (5.15) 

of the form (12.11 ) with s(x) = 1 + x 2 and b(x) = v + x for all x G I = R . In this case the diffusion 
in natural scale 

dX°{t) = (1 + (X°(t)) 2 ) dW°{t) , X°(0) = £ 

of (12.131 ) does satisfy F°(S° = 0) = 1 but does n ot satisfy the condition of (12.1 lb : it is a strict local 
martingale, studied for example in lCarr et al.l (120131) . We have 



„. . v + x . 1. . log(l + a^ 

H x =T— -2' Fx H^tan- 1 x +-^- 

1 + x z 2 



V{x) 



1 + 



in the notation of Subsection |3.1.1| Then ( 13.61 ) and Corollary 1 and Lemma 2 in lCarr et al.l (120 131) . which 
provide a distributional identity of the local martingale X°(-) in terms of the Brownian motion W°(-), 
applied with 



g{x) = cos(x + c)/cos(c) , f(x) = tan(x + c) , = (1 + x 2 ) 1 ^ 2 exp(z/tan 
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and C = 1, where c = tan yield the representation 



»(5 > T) = E° 
= E° 



1 + (X°(T)) 
l + £ 2 

exp (vW°{T) 



2\ 1/2 



^ exp fi/tan _1 pf°(T)) - i/tan _1 (£) 

J 1 {t°>T} 



T + u 2 T 



v 2 T 



with r° the first time that 14 / °(-) hits either 

o = a(f) := -(tt/2) -taiT 1 ^) or 6 = 6(0 := (tt/2) - tan" 1 ^). 

Alternatively, we observe that the equation (I5.15l l is again of the Lamperti-Doss-Sussmann form 
(15.171 ) with fj,(-) = v, and can thus be solved pathwise by the procedure outlined in Subsection 13.21 The 
end result 

X(t) = tan (W(t) + ut + tan" 1 ^)) , < t < S, 
identifies the explosion time as a first exit time for Brownian motion with drift 

S = m£{t > : W{t) + ut £ (a, b)} 



from the open interval (a, b). Formula (3.0.2) on page 309 of iBorodin and Salminenl (120021) now com- 
putes the distribution of this exit time as 

f°° { -v 2 t\ 

P(5 > T) = / exp ( exp(zya)e (i; b,b - a) + exp(i/6)6 (t; —a, b — a) )dt, (5.16) 



where 



<d)(t;u,v) := >^ ; — exp 



V2vrt 3 



(t; — u + 2fci>) 
2t 



denotes a certain inverse Laplace transform. 

In this example, the explosion time 5 has actually finite expectation. From Propositions 14 . 5 1 and 14 . 6 1 
the function on the right-hand side of (15.161 ) is the smallest nonnegative classical solution of the Cauchy 
problem 

dU , . (1 + x 2 ) 2 d 2 U , , , w 2 ,dU , . , . . 



with boundary condition W(0, x) = 1 for all x G 



□ 



Example 5.7. Diffusion with cubic variance function. In a similar manner, we consider the stochastic 
differential equation 



dX(t) = (X(t)) 



3/2 



dW{t)+(v+Ux(i)) 



1/2 



dt 



x(o) = s e i 



with state space I = (0, oo), for some given real number v. This equation is of the form (12.11 ) with 

s(x) = x 3//2 , b(x) = v + (3/4)x 1 / 2 , and the diffusion in natural scale 



dX°(t) = (X°(t)) 3/2 dW°(t), X°(0) = £ 
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of (12.131) satisfies the conditions of (12.10I ). We note that 



too 



v 3_ 

x 3/2 4 X ' 



F(x) 



-2v 



+ log(x 3 / 4 ), V{x) 



v 

T 



3.r 
32 



in the notation of Subsection l3-l.il so (13.61 ) gives the representation 



> T) = exp 



2v 

7T 



u 2 T 



E° 



X°(T) 



3/4 



cxp 



32 



X°{T)dt 



2u 



This probability can be computed explicitly, once the stochastic equation is cast in the Lamperti- 
Doss-Sussmann form (15.171 ) with /«(•) = v. Following the methodology of Subsection 13.21 we obtain 
the diffusion X(-) explicitly as 



X(t) 



o<t<s. 



This shows that the origin is inaccessible by the diffusion X(-), which can only explode to infinity, 
consistent with part (i) of Proposition 12.11 The explosion happens at the Brownian first passage time 



S = inf It > : W(t) + vt 



whose distribution is of course well known, namely 



P(5 > T) 



1/2 



ex P -7TT 



2* \Vt 



vt 



dt=: U(T,£); 



see also Section 3.5.C in Karatzas and Shrevel (1991). In particular, with v < we have P(5 < oo) = 
exp(4^/y / ^); whereas, with v > 0, we have F(S < oo) = 1. It is checked by direct computation, that 
the function U(-, ■) solves the parabolic partial differential equation 



dU 



(t,x) 



x 3 d 2 U 



(t,x) + \vx 



,3/2 



with boundary condition U(0, x) = 1 for all x € I. Indeed, from Proposition 
[/(•,•) is the smallest nonnegative solution of this linear parabolic equation. 



the above function 
□ 



A related example is discussed in Corollary 1 of lAndreasenl (120011) . 
Example 5.8. Bessel process of dimension 5 G (— oo, 2). Let us consider the stochastic equation 

5-1 



dX(t) 



2X(t) 



dt + dW(t), X(0)=£ 



with state space / = (0,oo) and dimension S G (— oo,2). The solution of this equation does not 
explode to infinity, but reaches the origin in finite time: P(5 < oo) = 1, where S = lim^oo S' n with 
S' n = inf{t > : X(t) < 1/n}. 

In our notation s(-) = 1 and, with v = 1 — <5/2>0, we have 



f(x) 



1/2 - v 



F(x) = log ( x 



A/2-v 



V(x) 



1/4 



2x" 
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in the notation of Subsection 13.1.11 The representation (13.61 ) then helps us compute the distribution of 
S as the expectation of a functional of the Brownian motion X°(t) = £ + W(t), < t < S° with 
S° = M{t > : £ + W(t) = 0}; to wit, 



P(5 > T) = E° 



Z°(T) 
X°(T) 



-2v 



-2v 



X°(T) 



u+1/2 



cxp 



1/4 



df 



-{S°>T} 



c, exp 



2T 



x 1 u exp 



2T 



dx 

, given in (15 - 8b - Here, 



(5.17) 
denotes again 



in terms of the modified Bessel function of the second type l v 
the probability measure under which the process X°(-) is a Bessel process of dimension 2v + 2 = 
4 — <5 > 2 ; the change of measure is proved in Exercise XI. 1.22, and the density of the process X°(-) is 
derived on page 446, of iRevuz and Yor dl999b . The above expression becomes exactly the same as that 
of (15.71 ) in Example l5.2l if one replaces the initial condition £ £ (0, oo) by its reciprocal. 

Now with the help of the monotone convergence theorem and of the substitutions z = x/ v2T and 
y = £ 2 / (2T), the expression in (I5.17t simplifies to 



s (5 > T) 



2y v/2 exp 
2y^ 2 



y) 

exp (-y) 

k=0 



z 1 v exp 



2 ) I v (2zjy) dz 



,v/2+k 



z l+2k exp 



-z 2 ) dz 



V 



v+k 



expH/) gr ( „ + t + i) 



1 



I» 



exp(-6)6 u - L d6 



Thus, we have 

P(5 > T) 



© < 



2T 



2T 



I» 



where 



e/(2T) 



exp (-0) 9 U - L d6 



H(u) 



u t l/ - 1 exp(-t) 



dt 



is the cumulative Gamma(^) probability distribution function, and the random variable (5 has Gamma 
distribution with parameter v . 



Of course, it is well known fr om the t i me-re versa l considerations i n Sect ion 2.1 of lGoing-Jaeschke and Yor 
d2003b - based on the results in lGetoor! (1979b and [pitman and Yon d 198 lb - that S has the distribution 
of £ 2 / (26); see also Section 14 in lKenl Jl982h . Here we just derived this fact from rather elementary 
Bessel process computations, without any need for time-reversal. 

It is veiy easy now to check "by hand" that the function (T, £) i— > F^(S > T) is a classical solution 
of the linear parabolic equation 



1 d 2 U , m ^ 5 



— -(T 

2£ oc ,ih 



(T,0 e (0,oo) x/, 



and satisfies also the initial condition U(0+, £) = 1 for £ G I and the lateral condition W(T, 0+) = 
for T G (0, oo). From Proposition 14.61 it is the smallest nonnegative classical solution of this initial / 
boundary value problem. □ 
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Remark 5.1. A generalization of Example \5.8\ Barndorff -Nielsen et al. construct, for a given set 

of parameters, a one-dimensional diffusion X(-) on the interval / = (0, oo) such that the corresponding 
time to explosion S has a generalized Gamma distribution. If the parameters are chosen so that the 
generalized Gamma distribution is exactly a Gamma distribution, then their construction yields the Bessel 
process of Example 15.81 □ 



Example 5.9. Explosion to infinity. Let us now consider the situation with / = M, s(-) = 1 and 
b(x) = exp(/3x) for some f3 > 0. In this case the Brownian motion X°(-) = £ + W°(-) has no 
explosions, i.e., P°(S° = oo) = 1; the process X(-) with dynamics 



(5.18) 



X(-) = £ + / exp (pX(t)) dt + W(-) 
Jo 

explodes (to +oo) in finite time, that is P(S < oo) = 1; and 



»(5 > T) = E° 



,2«A-(<) dt 



holds for all T G (0, oo). From Formula (1.30.7) on page 196 in iBorodin and Salminenl (|2002h . we 
obtain the distribution of the explosion time as 



P(5 > T) = /3 2 



oo roo 



oo JO 



exp ((e"* - J*) ft - (Py/2)) 



exp 



8sinh(/3y/2) V Pt&nh(Py/2) J 

^ {p S ,n H Py/2)) dydZ > 



where 



i t (z) := CJ 1 (i^z)) 



z exp(vr 2 / '(At)) 



exp I — zcosh(u) — — ) sinh(it) sin ( — ) du 
o V At J V 2t 



denotes an inverse Laplace transform related to the modified Bessel function of the second type in ( I5.8I ). 
The substitutions of Py/2 by ijj, and of exp(/3z/2)/y / 7J by (, and the change of variable 



exp (Pj/2) 
VP 



transform this expression as 

2exp(2vr 2 /(/3 2 T))v / 2r ? 



P(5 > T) 



ttPVttT 



OO /"OO /"OO 



exp (C 2 — rj 2 - i[)) 
o ./o ./o (sinh(V')) 2 



tanh(^) J 



( > cosh(n) 2m 2 \ . . /4ttu\ 



du dip dC ■ 



□ 
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